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QJ , Abstract 

, The theory of quadrature domains for harmonic functions and the Hele-Shaw problem of the 

' fluid dynamics are related subjects of the complex variables and mathematical physics. We 

I present results generalizing the above subjects for elliptic PDEs with variable coefhcients, 

■ emerging in a class of the free-boundary problems for viscous flows in non-homogeneous 
^ I media. Such flows posses an infinite number of conservation laws, whose special cases 

may be viewed as quadrature identities for solutions of variable-coefficient elliptic PDEs. 
I If such PDEs are gauge equivalent to the Laplace equation (gauge-trivial case), a time- 

■ dependent conformal map technique, employed for description of the quadrature domains, 
I leads to differential equations, known as "string" constraints in the theory of integrable 

|0 ■ systems. Although analogs of the string constraints have non-local forms for gauge-non- 

I trivial equations, it is still possible to construct the quadrature domains explicitly, if the 

• elliptic operator belongs to a class of the Calogero-Moser Hamiltonians. 

■ 1 Introduction 

^ . Recently, the constant-coefficient Helc-Shaw problem has received a good deal of attention due to its 

^ ' connection with the theory of integrable hierarchies in the dispersionless limit [8], [15]. These integrable 

• i-H . structures turned out to have natural interpretations in the classical free-boundary problems. 

' Although concerned with various aspects of integrability, the present paper is primarily devoted 

^ . to construction of solutions for a more general class of the variable-coefficient free-boundary prob- 

■ " " ' lems. Namely, motivated by practical applications, we derive a class of rational solutions of the 

variable-coefficient Hele-Shaw problems, related to the dihedral Calogero-Moser systems. These solu- 
tions provide examples of quadrature identities for elliptic PDEs with variable coefficients. The theory 
of quadrature domains for such PDEs is highly reminiscent of that for harmonic functions. 

The structure of the paper is as follows: A brief summary of the theory of the Hele-Shaw flows, as 
well as that of the quadrature domains, is given in the next three sections, followed by introduction 
of the time-dependent conformal map technique in Section 5. We digress into consideration of the 
gauge-trivial problems, considered in the above mentioned works on relationship between dispersionless 
integrable hierarchies and the Hele-Shaw problem, in Section 6. Sections 7 and 8 are devoted to 
derivation of explicit results for a class of the gauge- non-trivial problems. Concluding remarks are 
given in Section 9. 
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2 Variable Coefficient Hele-Shaw Problem 



The Hele-Shaw flow governs the dynamics of the boundary = in the plane separating two 

disjoint, open regions Q = Q{t) and C\fi in which scalar fields are defined. These may be interpreted as 
the pressure fields for two incompressible (viscous and non- viscous) immiscible fluids, trapped in a gap 
between two plates or propagating in a layer of the porous medium. If the size of the gap or the layer 
height is negligible, the velocity of the viscous liquid can be averaged across the direction perpendicular 
to the layer or the plates respectively, and the problem becomes effectively two dimensional. The 
velocity of the viscous liquid is then proportional to the two dimensional pressure gradient (the Darcy 
law) 

V = -VP (1) 
The flow is incompressible and the liquid velocity satisfies a continuity equation 

V • u = (2) 

The region can be chosen to be bounded and occupied with a viscous liquid, so that (1) holds in that 
region. It will be referred as the "interior" region. It is surrounded by a non-viscous liquid occupying 
the "exterior" (unbounded) region C\^l. Since the flow of the liquid occupying the exterior domain is 
inviscid, it must be driven by the pressure field with vanishing gradient and, therefore, the pressure 
is constant in that region. Without loss of generality the pressure can be set to zero in C\il. It is a 
continuous function across the moving boundary dQ{t), so that 

P(dn) = (3) 

and the normal velocity of the boundary coincides with that of the viscous flow at dVl 

n ■ ^dil = n - V (4) 
dt ^ ' 

where n denotes the outward normal to the boundary. In the present and the next section we consider 
flows that are driven by a single point source located inside O. Without loss of generality we can locate 
it at the origin z = ^ 

P — > — log(2:^), as z ^ 

where z := X + iY, and X, Y stand for the Cartesian coordinates on the plane. The last equation, 
together with the Darcy law (1) and the continuity equation (2) leads to 

AP = -'k5{X)S{Y) 

where 6 denotes the Dirac delta-function. In this paper we deal with a generalization of the above 
problem, namely the Hele-Shaw problem with coefficients depending spatial variables X and Y 

V = -kVP, V • (riv) = (5) 

where k = K(X,Y),r] = rj(X,Y) are arbitrary functions of X,Y, sufficiently regular in O. As in the 
constant coefficient case, the pressure has a logarithmic singularity at the origin, i.e. 



when the flow is driven by a single point source. 

Note, that one has to modify the asymptotic condition (6) if kt] vanishes or has a singularity at 
the point source, or if several point sources coalesce in a special way ("multipolc" sources). 

We consider a situation when the exterior region is occupied by a non-viscous liquid, so that (3) 
holds. Equations (3), (4), (5), as well as asymptotic conditions (e.g. (6), for a single point source) set 
the Helc-Shaw problem with variable coefRcients, describing, for instance, propagation of the liquid in 
a thin horizontal layer of non-homogeneous porous medium, whose permeability and height depend on 
X,Y. 



3 Conservation Laws 

Prom (5), (6) it follows that the pressure satisfies the second order elliptic PDE 

(Vk77V)P = -Tr6{X)S{Y) (7) 

that replaces the Laplace equation of the constant coefflcient problem. Let (p{z, z) be any regular in 
the interior domain solution of the equation 

(V.,V)0=(|«4 + |.4)^ = O (8) 

Introduce the following quantities 

M[(l)]= [ r](t)dXdY (9) 

Jn(L) 

and estimate their time derivatives. Considering an infinitesimal variation of the interior domain 
Cl{t) — > Cl{t + dt) along the boundary dQ,{t), we get 



dM[(i)] 
dt 



rjcpVndl 
dn{t) 



where dl is the boundary length element and w„ is the normal velocity of the boundary. From (4) and 
(5) it follows that Vn = v ■ n = —Kn ■ VP. Since P{dQ,) = 



dM[4>] 



{PKrjVcj) ~ (l>nriVP) ■ ndl 

dQ{t) 



dt 

Applying the Stokes theorem and taking (7) and (8) into account, from the last equation, we get 

dM[(l)] 



dt 

that, according to (6), becomes 



(PKrjVcj) — (j)K.r]\/P) ■ ndl 

z=0 



^^[^]=.0(O) (10) 



dt 

Therefore, the quantity M[(()] is conserved for any (regular in O) solution of (8), such that (p{0) = 
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The Richardson harmonic moments /q^j-j z'^dXdY [10] of the constant coefBcient Hele-Shaw prob- 
lem correspond to the special case 



K^rj=l, (f){z,z) = z'', k>0 

To the author knowledge, the conservation laws for the variable-coefficient problem were first presented 
in [3]. 

4 Quadrature Domains 

The derivation of conservation laws can be easily generalized to the case of several sources. Consider 
the flow driven by N sources of the time-dependent powers qk{t) that are located at points Zk,k= 1..N 
in fl 

— Qk{t) 

P — -— -log|z-Zfe|, as z^Zk, Zk = Xk+iyk&^ 

Then by arguments, similar to those used for the single point source problem we arrive at 

N 



dM[(j)] ^ , , , , _ , 
=Tr ^qk{t)(p{zk,zk) 

fe=i 

for any (j){z,z), regular in O and satisfying (8). It follows that 

M[(P]{t) = M[(t)]{0) +Try2Qk<l>izk,Zk), Qk = Qk{t)= / qk{t')dt' 

k=i Jo 

From the last equation, we see that M[(j)\{t) does not depend on history of sources and is a function 
of total charges Qi, produced by the time t. This fact reflects the integrability of the problem, where 
flows produced by different sources commute. 

Consider the special case when M[(/)](0) = 0. It describes the injection of the fluid to an initially 
empty Hele-Shaw cell. Taking into account the definition (9) of M[(/)] we obtain 

ri(j>{z, z)dXdY = Try2 QkHzk,Zk) (11) 
-^^ k=i 

This is an identity expressing the integral over the domain Q in the left-hand side as a sum of terms 
evaluated at a finite number of points, given on the right-hand side. The special case rj = k = 
1 provides quadrature identities for harmonic functions. Special domains fi, possessing the above 
property, are called quadrature domains [11], [12], [14]. The simplest example of a quadrature domain 
is a circular disc, produced by a single point source in the constant coefficient Hele-Shaw problem. 
The corresponding quadrature identity is a "mean value" theorem for harmonic functions. Equation 
(11) is a generalization of quadrature identities appearing in the theory of harmonic functions to the 
case of elliptic equations with variable coefficients. The quadrature domains for such PDEs are, thus, 
solutions to the variable coefficient, interior Hele-Shaw problems with zero initial conditions. In the 
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sequel we will mainly deal with situation when groups of sources coalesce in such a way that (11) 
becomes 

JV 

/ r]^{z,z)dXdY = 7ry2Qk[(t>]izk,Zk) (12) 
fc=i 

Where Qk is a finite-order, differential in z, z operator of the following form 

0. = qt + E te^li + Q^%) , = or, (13) 

For such combinations of multipolc sources, the integral over in the left-hand side of (12) is expressed 
as a sum of terms involving values of function (f) as well as a finite number of its derivatives at a finite 
number of points inside O. The boundary dQ{t) of the quadrature domain is a solution to the variable 
coefficient Hele-Shaw problem with pressure satisfying the following equation 

fc=l i=l 

Note, that operators Qk do not contain mixed derivatives gfrrgp^, since, due to (8), is a linear 
combination of the first order derivatives of (j). The coefficient in front of ^ must be complex conjugate 
of that in front of since both (f>{z, z) and its complex conjugate satisfy (8). 



d' 



at dzi 



dt dzi 



5 Time-Dependent Conformal Maps 

Usually, a time dependent conformal map technique is implemented to find explicit solutions to the 
constant coefficient Hclc-Shaw problem and, in the special case of zero initial conditions, it is also 
an efficient method of constructing quadrature domains. This technique is also applicable to the the 
variable coefficient case, leading to explicit solutions of some non-trivial problems. 

Introduce a "mathematical" u>-plane and denote by z{w,t) a conformal map from the unit disc 
\w\ < 1 in the w-plane to a simply-connected interior region ^l{t) in the physical z plane. According 
to the Riemann mapping theorem a one- to one analytic in \w\ < 1 map 

z{w,t) =r{t)w + ^u^{t)w^+'^ (14) 

j>0 

exists, such that the unit circle in the w-plane is mapped to the (analytic) boundary contour in the 
z-plane 

z{w,t) edQ.{t), if \w\ = l (15) 

Alternatively to the derivation of Section 2, one can estimate the time derivatives of M[<f)], transforming 
the two-dimensional integrals (9), taken over O, to line integrals along the unit circle in the t«-plane. 
Introducing a function ^(z,z) such that 

n{X,Y)(l>{z,z) = —az,z) (16) 
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by the Green theorem and (15), we may rewrite (9) as 



M\(j>] = ^ i ^dz = ^ i i^dw (17) 
Note that r(t) in (14) can be made to be real and 

w = llw, -z=-z[l/n^,t) = '-^+Y^^ (18) 

i>0 

along the boundary. In the sequel wc mainly deal with z{w,t) evaluated at the boundary w = l/w 
and therefore we use z to denote z{l/w,t) (if not otherwise specified). It follows from (16) and (17) 
that 

dM[0] _ 1 / ifdzdz dz dz\ 1 / d ( dz" 



dt 2i /|u)|=i ydw dt dt dw ) 2i J\w\=i dw y dt _ 

which equals 

dM[^] _ 1 r , {z{w,t),z{l/w,t)} 



dt 2i 



w\ — l 



W 



provided ^ in (16) is univalent at dQ. In (19), {, } denotes the Poisson bracket 

f I -.gf (.0, 

defined on the cylinder parameterized by coordinates {w = e^^,t) with real 6 and t. By virtue of (19) 
the evolution equation (10) may be rewritten as 

J- I ^r,i^dw = 0(0) (21) 

6 Gauge Trivial Problems, String Equations 

Consider equation (8). Its solution (/> also satisfies 

m=0, L:=i^-^—^^—+^-^—^^ — , ^z.2) = i,{-z,z) (22) 

The elliptic second order differential operator L is amenable, by a gauge transformation, to the two- 
dimensional zero-magnetic field Schroedinger operator 

H:=-i;L^ = ~^ + V{X,Y), V=^, //[V<^] = 0, if L[</.] = (23) 

The potential V vanishes when tp = h, where h stands for a harmonic function 



L = £o := h~^Ah, £o 



= 0, A;>0, h:=^{z) + '^{z) (24) 
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The elliptic operator Co is gauge equivalent to the Laplace operators. We call such operators and 
corresponding Hele-Shaw problems "gauge- trivial" . According to (9), (24) the following quantities 



= / ^dXdY, Ik= I ^z^dXdY, k>0 
Jn(t) h jQit) h 



/Q(t) JQ{t) 

are linear and constant in time. 

Using (24), from (21) we see, that in the gauge-trivial case 



f Tj f Z^ T) 

(p — — {z,z}dw= d) — — {z,z}dw = 2niSko, k>0 
Jm=i w h w h 



— — {2, z}dw = f 
i\w\=i w h 

or, taking (14), (18) into account 



/ iT{z,z} -A—dw = i (^{z,z} -l]—dw = Q 
J\w\=i J w J\w\=i ' w 

The expression in the parentheses of the last equation must vanish, since jw = r^w^^^ + ... and 
^fc ^ ^k^^j-k-i _ fc > form a basis of an arbitrary Laurent series in w. Therefore, z, z satisfy the 
following differential equation 

{z{w,t),z{l/w,t)} = h/r] 

or equivalently 

{q{z,z),q{z,z)} = l, (d,q){d,q) - {OMd-zq) = h/r^ (25) 

known as the "string" constraint in the theory of the dispersionless integrable hierarchies [5], [8], [13], 
[15]. In the special case /i = 1,7? = 1, or equivalently q{z,z) = z, corresponding to the constant 
coefficient Hele-Shaw problem, the conformal maps satisfy the following constraint 

{z{w,t),z{l/w,t)} = l 

known as a Galin-Polubarinova equation [4] , [9] in the theory of the Hele-Shaw flows. 
The string equation (25) is preserved by the Lax-Hamilton flows 

fj 0' fly 

^ = {z,'H}, ^ = {z,'H}, H = H{w,t,T) 



Indeed, 



^{q{z, z), q{z, z)} = ( ^^i^, ^-(z, z)] + \q{z, z), \ = {H, {q{z, z), q{z, z)}} = 



dr [ dr J [ dr 

i.e. functions z{w, t), z(l/w, t), satisfying (25) belong to invariant, under the action of the Lax-Hamilton 
vector fields •}, subspace of space of functions of w, t. It is also necessary to satisfy the condition of 
the form-invariance oi z{w, t) along the Lax-Hamilton flow lines, i.e. such Lax-Hamilton functions must 
be chosen, that z{'w,t) will remain the Taylor series (14) along the corresponding flows. So selected 
Lax-Hamilton functions, generate symmetry transformations of the Hele-Shaw problem, mapping con- 
tinuously one solution of the problem into the others. An ableian subset of these transformations. 
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forms a parametrizable set of deformations, leaving ivariant equations governing the Hele-Shaw flow. 
The Lax-Hamilton functions generating such abelian subset can be conveniently chosen as 

Hk^{z-'')<o + l/2{z-%, iffc = (z-'=)>o + l/2(z-'=)o, A:>0 

where ()>oj ()<0) ()o stand for negative, positive and zero parts of the Laurent expansions 

(/)>o 

fe>0 

The flow equations 

dz 



constitute the two-dimensional Toda hierarchy in the dispcrsionless limit (2dToda) or Sdiff(2) hierarchy 
[13]. Due to commutativity of the 2dToda vector fields, the maps (14), (18) are functions of the 
deformation parameters ti,T2, .... fi,f2, ... (the 2dToda "times"). The 2dToda system is an integrable 
Hamiltonian system of PDEs for the coefficients r(t, ti, fi...), Ufe(t, n, fi...), f2fc(t, n, fi...), fc > 
of the scries (14), (18), obtained by equating (26) as Laurent series in the dummy variable w. The 
string equation (25) defines a reduction of the 2dToda hierarchy. The connection between the 2dToda 
hierarchy and the Hele-Shaw problem was first found in [8], [15]. 

In the context of the constant-coefficient Hele-Shaw problem, the 2dToda "times" r^, can be nat- 
urally interpreted as harmonic moments of the domain fi. More precisely, the fcth harmonic moment 
Ici(t T f ) z'^dXdY evolves linearly in rj, and is constant along the other 2dToda flows [5], [8], [15]. 



(/)o:=/o, (/)<o:=E^"^'' if / = E/' 



k<0 



fcez 



{nk,z} 
{nk,z} 



dz 



^{nk,z} 



k>0 



(26) 



7 Gauge Non-trivial Problems, Quantum Integrable Systems 
on Plane 

In the gauge trivial cases, equation (21) can be transformed into the differential "string" equation (25) 
for time dependent conformal maps. In contrast to the gauge-trivial case, similar differential represen- 
tations of the gauge-non-trivial integral relations (21) seem to be generically impossible. Nevertheless, 
we can still construct explicitly the quadrature domains for special gauge-non-trivial elliptic PDEs. In 
this section, we derive sets of solutions to (8) and related conserved quantities that will be used for 
construction of such quadrature domains. 

An explicit evaluation of conserved quantities is possible when a gauge-non-trivial elliptic operator 
L in (22) is equal, up to a gauge transformation, (or (23) is equal) to a Hamiltonian of a quantum 
integrable system on the plane. 

We consider problems (5) with such kt] that the corresponding second-order elliptic differential 
operator L (22) can be related to the Laplace operator by a differential operator T 

TA = LT (27) 

T is usually called the intertwining operator and (27) is an intertwining identity. The simplest in- 
tertwining identity corresponds to the gauge trivial case (22) when T = Tq := is a zero-order 
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differential operator and L = Cq := j^^h. If, fiowever, T is a differential operator of a non-zero order, 
the corresponding L equals, up to a gauge transformation, a Hamiltonian of a non-trivial integrable 
quantum system on the plane. Indeed, from (27) it follows that any eigenfuction '^{z, z, A) of A 

A[^{z,z,X)] = X^{z,z,X) 

is transformed by the action of T, '^{z, z, A) — > T[^{z, z, A)], to an eigenfunction of L having the same 
eigenvalue A or to zero. 

We start with a class of examples, in which k, t] vary only in one direction (say X-direction) and 
equal 

'^=Y^' ^ = 1' " = 0'1.2,- (28) 
The corresponding elliptic operator (22) 

did 

L = Cn:=X^^ — ^— + g^ (29) 

is amenable, by a gauge transformation, to the two-dimensional Schrodinger operator 

d cP" ti(ti -\- 1) 

Hn = = -gy2 + Qjl2 ^^^^ 

which , when rewritten down in the polar coordinates, has the following form 

The operator 5„ admits the following alternative factorisations 

Sn = {-§Q + ntan(6')) (^ - ntan(6')) - = 

(^ - (n + 1) tan(e)) (^ + (n + 1) tan(^)) - (n + 1)^ 

and therefore 

(n + 1) tan(^)^ 5„ = 5„+i + (n + 1) tan(e)^ 
In view of (29), (31), this leads to the intertwining identity 

f d d Y\ f d d Y\ 

T„A = r„T„, T„ = X- (X- -Y^Hn + D^j ' ' ' (^^ " + ^) (32) 

Note, that the intertwining operator is not unique. To see this on the example of equation (29), we 
can use the translational invariance of the latter along the F-direction, shifting Y in T„ by a constant 
A, and so obtaining a linear combination of operators 



d 



T„ + ^A^r, 



(i) 

'\ J. 5 
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Tn as well as T„, are all intertwining operators of the nth order. They are homogenous polynomials 
in. X,Y,-^,^. For instance, in the simplest n = 1 case 

T.=X[X^-Y^)-Y, T«=^^-l (33) 

It is easy to see that operators T^"' (e.g. t'^^ in the above example) are X-dependent only, i.e. 

ri") = ^"^ + E«'=;"^'a^ (34) 

fe=0 

Operator (34) can bo also obtained alternatively by exploiting separation of (29) in Cartezian coordi- 
nates X,Y, through a chain of factorizations leading from to -^^" ^ jfW^- 

Images of all intertwining operators, acting on harmonic functions, coincide for a given n. Therefore, 
one can choose any operator from T„, Tn\i = l..n (or their linear combination) to construct same set 
of solutions to (29). In so doing, we return to the n = 1 example (33), choosing t'^^ In this example, 
it is convenient to make a shift of z by the distance x\ along the X direction, displacing the singular 
line X = of the problem (29) to X = —x\. Then, without loss of generality we can locate a source 
at point z = Q. It follows that the functions 

<^fe;i = 2T[%''+^] = {k + l){z + z + 2xi) z'' - 2z''+\ fc = 0, 1, 2, ... 
and their complex conjugates, form a set of solutions of equation 

^'' + '^^ dxJxT^dx + W^) ['^^^^] = ° 

related to the variable-coefficient Hele-Shaw problem with k = 1/{X + xi)^, rj = 1. Plugging this 
set into (9), by (10) we get an infinite number of quantities 

Mk-i = M[cPk;il (35) 

that are linear and constant in time. They form a complete set of local coordinates. A simple way to 
see the latter is to tend xi to infinity 

- z'^ + '— > z", as a;i ^ oo 



2(fc + l)xi 2(fc + l)xi 

observing that the set (j)k;i tends continuously to a basis z'',k = 0,1,.. of functions analytic in a 
neighborhood of z = 0. The linearizing coordinates (35) are, therefore in one to one correspondence 
with harmonic moments J^^ z^dXdY , at least in some neighborhood of infinity. The corresponding 
variable-coefficient Hele-Shaw problem also transforms continuously into the constant coefficient one. 
Since harmonic moments are local coordinates for a generic set of simply-connected domains O (see 
e.g [14] and references therein), so are M[(j)k], fc = 0, 1, .... 
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A similar argument about completeness may be applied to the rest of examples considered below. 
For instance, it is convenient to use the following set of regular solutions to (29) when dealing with an 
arbitrary n problem for a flow, driven by a finite combination of multipole sources located a.t z = zi 

.^^") = i, <p'c\z,z)=T(r^[{z-z,r+% ^'c\z,z)=T(r\{-z--z,r+% k>i (36) 

The following class of variable coeflacient Hele-Shaw problem, whose conserved quantities write in 
terms of polynomials in z, z, generalizes (28) 



1 



K = TTT^ rrrr, ?? = 1, n>Z>0 (37) 



An elliptic operator (8), corresponding to (37) 

r-r o ^""^ /^-m-i 9 d \ 2lm / d d \ 

5^5^ ^- + ^-V dz^ dz) z--z™V 9z dz) ^ ' 

equals, up to a gauge transformation, the Schrodinger operator of the Calogcro-Moser system related 
to (dihedral) group of symmetries of a regular 4m-polygon (2m-polygon if Z = 0). 

nn,l;m = {z"" + Z^Y {z"" - zT)' lin,l;m{z^ + 0'")-"(z'" - z'")"' = ^ + f ^ + 

Q _ _ m^n(ra+l) _ m^l{l+\) 

On,l;m — -gff! cos(m6>)^ sin(m6i)^ 

By analogy with Sn in (31), Sn,i;m admits factorizations leading to the intertwining identity 
where the intertwining operator Tn,i-m can be represented in a form of the Wronskian [1] 



-l-n,l;m[T\ = P , Q\r,(r,-U/2 ■ ( /l> "> /fc •= ^^t 

cos(m6')"v" ^)/^ sm(m&)H' 



(39) 



l<i,j<k 



with 

sin(fc(m6' + 7r/2)), k = l,2,..n-l 

cos {{2k + I - n){m0 + 7r/2)) , k = n - I + l,n - I + 2, ..n 



Tn,i;m is a differential operator of the nth order and is a homogenous polynomial in z, z, It 
transforms any holomorphic function f{z) into a polynomial in z, which is crucial for construction of 
quadrature domains. 

It is noteworthy that the intertwining operator Tn^i-m is also not unique. Although, in difference 
from the I = 0,m = 1 case we cannot apply a simple argument connected with the translational 
invariance of the system along the K-direction, another intertwining operator, constructed by the 
Dunkl method [2], exists. Various operators connecting A with Cn.i-.m, that are compositions of the 
intermediate intertwining operators constructed by the Dunkl method and those constructed through 
factorizations, can be also obtained. Similarly to the Z = 0,m = 1 case, images of all so obtained 
intertwining operators, acting on harmonic functions, coincide for fixed n,l, and to derive a complete 
set of solutions to (38) one may use just one of them, say (39). 
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8 Construction of Gauge-non-trivial Quadrature Domains 

Let us start from gauge-non-trivial examples of the Hele-Shaw flows (5), (28), that are driven by a 
multipole source located a.t z = zi. Since the quadrature domain fl can be viewed as a solution to 
the Hele-Shaw problem that develops continuously in time, starting from zero initial data, its form is 
defined by the condition that the quadrature identity (12) holds for any solution of (29), regular in O. 
We show first, that curves, parametrized by polynomial maps of any non- negative degree s > 

s 

z{w) = zi+rw + '^Uk-iw'^, z{e''^)edn, (40) 

fe=2 

where r can be made to be real, are boundaries of the quadrature domains for solutions of (29). As 
in the constant-coefficient problem, forms, sizes and positions of these domains are functions of s -|- 1 
free complex parameters zi,r,Uk,k = 2..s. 

As follows from the previous section, any solution to (29), that is regular in O can be represented 
as (34) 

^ = 2T(")[/(.)] = f:^(. + .-)^^, a„;„ = l 

k=0 

where f{z) is analytic in vicinity oi z = z\. Substituting this solution into (12) and using the Green 
theorem, we can transform the integral over O in the left-hand side of (12 ) to the line integral 



Using (40) and taking analyticity of f{z) into account, we get 

(n+l)(s+l)-2 ^ 

d 



I 

Jn 



(j)dXdY = ^ Ui{zi,r,ui,..,Us,ui,..,Us)i ■ ) (41) 



i=0 



Comparing (41) with the right-hand side of (12), we see that the latter must contain s{n + 1) — 1 
derivatives of and is equal to 

(s+l)(n+l)-2 . 



j=0 



where, Vi are linear functions of \ Q^°\ Eqs. (41) and (42) lead to the 

non- homogenous over-determined linear system of 2(s -|- l)(n -|- 1) — 1 equations 

Ui = Vi, Ui = Vi, i = 0..(s + l)(n + l)-2, gf) = 0l°) (43) 

for 2s(n + 1) unknowns Q'^\ k ~ 0..s(n + 1) — 1. Although the number of equations exceeds the 
number of unknowns by 2n+ 1, not all these equations are independent. For instance, solving (43) for 
circular domains, s — 1, z = zi + rw, we get the simplest gauge-non-trivial quadrature identities 



4 

7rr 



(j)dXdY = nr (j){zi) + - — ^ + ^ , zi=xi+iyi 
\z-z,\<r 4a;i \dz dz J 
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for solutions of (29) with n = 1 and n = 2 respectively. Moreover, a circular solution to the arbitrary 

n variable-coefficient Hclc-Shaw problem (28) exists, if the flow is driven by a special combination of 
a monopole, dipole ... n + 1-pole sources located at the same point zi. 

To show that the polynomial solutions of an arbitrary non- negative degree exist for any n in (28), 
(i.e that the number of independent equations in arbitrary n, s case equals the number of unknowns 

Qi\Qi'' and these equations are compatible), it is more convenient to check the quadrature identity 
(12) on a complete set (36) of solutions to (29). 

According to (41), (42) and (36), we have to evaluate 

(s+l)(n+l)-2 (s+l)(„+l)-2 

3=0 i=o 

on f{z) = {z — 21)"+*^ (and f{z) — {z — £1)"+*^ respectively) at point z — zi for fc > 1. Also there 
is one more equation, obtained by substituting = = 1 in (12). The latter is equation for Q^i\ 
that, as seen from (12), has a real solution, and equation = Qi'^ in (43) is satisfied. So, there 
remains 2s(n + 1) — 2 unknowns Q'f \ Q'i\i — l-.s(n + 1) — 1. 

The highest derivatives in (44) are of the order (s + l)(n + 1) — 2, so that {z — ZiY (or {z — zxY 
for the second equation in (44)) is annihilated at ^ = 2:1 if j > [s + l)(n + 1) — 2. As a consequence, 
(44) is identically satisfied for k > s{n + 1) — 1 in (36), and there remains 2s{n + 1) — 2 equations and 
the equal number of imknowns Q'f KQ'f ^.j = l..s(n + 1) — 1. 

We now have to prove the compatibility of remaining equations. (44) is a non- homogenous system of 
linear equations for unknowns Q^i\Q^i\Q^i \ fixing their dependence on parameters Zi,zi,r, Ui, Ug, 
ui, Ug. Such a system is compatible if its determinant docs not vanish i.e. its homogenous part docs 
not have nontrivial solutions. Let us suppose that it does. Remind that the homogenous part of the 
system has been obtained by the action of the operator 



s{n+l)-l , ,j „j 



in the right hand-side of (12), to an arbitrary solution of (29) at point z = z\. So, if the determinant of 
the system vanished, then there would exist such , Q^^^ , Q^^^ , A; = l..(s + l)n — 1, that the operator 
(45) would annihilate any solution of (29) at 2; = ^i, i.e. 

Qiri")[/(z)] =0, at z = zx (46) 

where /(z) is any analytic in 17 function. If the above were true, then changing continuously the 
position z = zxoi the source, we could pick such coefficients Q^^^ = Qf\z, z), Q^^ = Qi\z, z), Q^^ = 
Qi\z,z), that (46) would hold at each point in some region of the plane. In so doing we could 
construct a such operator Qi, with coefficients depending on z, that QiTn"'\f{z)] = in some region 
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of the plane for an arbitrary f{z). But this is evidently impossible, since the highest symbols of Qi 
and T^"^ contain pure derivatives in z, so does their composition. 

Thus, the system of equations (43) is compatible and has a unique solution 

Q^'-' = Qf\r,ui, ..,Us,ui, ..Us), Q^i^ = Qf\xi,r,u\, ..,Us,ui, ..Us), k = l..s{n+ 1) - 1 

defining quadrature domains with boundaries parametrized by polynomial conformal maps (40) of 
an arbitrary non-negative degree s > 0. These domains are solutions to the Hele-Shaw problems, 
describing the free-boundary flows driven by a multipole (combination of a monopole, dipole, 
s{n + l)-pole) source located at point z = z\. 

Analogous analysis can be applied to the problems related to general dihedral systems (37), (38). 
It also leads to the conclusion about existence of polynomial solutions (40) that depend on s + 1 free 
parameters. The flow, in this case, is driven by a combination of a monopole, dipole, ..s{m(n + /) + 1)- 
pole sources located at z = zi. 

9 Conclusion 

We have shown that, similarly to the constant coefficient Hele-Shaw problem, a class of the variable 

coefficient problems, connected with the dihedral quantum Calogero-Moscr systems, admits polynomial 
solutions for the flows driven by a finite number of multipole sources. In particular, the spaces of 
polynomial solutions have same dimensions in both the constant and variable coefficient cases. 
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